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ABSTRACT 

By generalizing earlier work of Johnston & Kulkarni, we present a detailed descrip- 
tion of the reduction in the signal-to-noise ratio for observations of binary pulsars. 
We present analytical expressions, and provide software, to calculate the sensitivity 
reduction for orbits of arbitrary eccentricity. We find that this reduction can be quite 
significant, especially in the case of a massive companion like another neutron star or a 
black hole. On the other hand, the reduction is less for highly eccentric orbits. We also 
demonstrate that this loss of sensitivity can be recovered by employing "acceleration 
search" or "acceleration-jerk search" algorithms. 
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1 INTRODUCTION 

It is now well established that studi es of binary pul sar systems provide a unique means to test relativistic theories of gravity in 
the strong-field regime (|Stairaii2005i : iKramer fc Stairs»,2008, l. with the most stringent tests often possible in systems with short 
orbital periods and high eccentricities. Detecting such systems in pulsar searches is challenging due to the Doppler modulation 
of the pulsed signal as the pulsar moves with respect to the centre of mass of the binary system during an observation. In 
an earlier study, Johnston & Kulkarni (1991; hereafter JK91) derived an analytical framework to compute the reduction of 
the signal-to-noise ratio due to binary motion for the case of circular orbits. This work also considered the improvement in 
sensitivity due to "acceleration searches" in which the data are corrected for binary motion either in the time ( Andersonlll993l ) 
or frequency domain ( Ransom, 200 1). 

Although circular orbits are applicable to many known binary pulsar systems (Lorimei 20081 ). for the purposes of a 
growing number of studies, it is desirable to compute the loss of sensitivity in eccentric binary systems for different types 
of companions. The purpose of this paper is to extend the work of JK91 to consider the case of orbits with arbitrary 
eccentricity. This framework allows us to compute the sensitivity degradation for any type of binary system, therefore allowing 
detail ed investigation of survey sen sitivity limits, which is app licable to a number of statistical studies (see iRidlev fc Lorimei 
I 2OI0I ) for example). While Ramachandran fc Portegies ZwartI l| 19981 investigated this, they did not provide explicit analytical 
expressions, nor did they extensively explore the parameter space implied by pulsar surveys and the different types of binary 
systems. We will discuss more about their work later when appropriate. 

The plan for the rest of this paper is as follows. The analytical expressions to describe this problem are developed in 
Section [2] numerical methods to employ those expressions are elaborated in Section [S] and some applications of the model in 
different cases are reported in Section [H Finally, our conclusions are given in Section [S] 
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2 FORMULATION 

Following JK91, we consider a pulsar with spin period Pp. The signal emitted from the pulsar at an arbitrary time t can be 
written as a summation of Fourier components: 

fc — oo 

'S'p(i) = ^ akex.p{ikijjpt + i^k) , (1) 



where ujp = 2-k jP^ is the angular spin frequency of the pulsar, tph is the phase factor of the A:*'' component and i 
signal received on the Earth 

k — oo f / , 

a 



Snii) = ^ «fc exp ( jfctjp ( * + ~ ) + 



-1. The 



(2) 



where d is the distance between the pulsar and the Earth at time t and c is the speed of light. The m' ^ harmonic of the 
received signal is 



Sii{t) = a,„ exp I imujp lt-\ — ] + itp' 



Using a Taylor series expansion of d, we can write the distance to the pulsar 



d = do + viot + 



aiot jiot 



+ 



(3) 

(4) 
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where vi is the line-of-sight velocity, ai is the line-of-sight acceleration, and ji is the line-of-sight jerk at the time t; vio, aig 
and jio are the values of these quantities at t = 0. do is the distance between the pulsar and the Earth a,t t — 0. An isolated 
pulsar moves with a constant velocity, so a;, ji (and higher order terms) are zero; but, this is not the case for a binary pulsar. 
Similarly, for vi, the Taylor series expansion gives 

j-2 

(5) 



vi = vio + aiot + + , 



Integrating this expression with respect to time, we get 
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So Equation Q becomes 



d = do + vi dt 



(6) 
(7) 



To detect a pulsar, one performs a Fourier transform of the time series of the received signal as T [SR{t)] — j]^ S'jj(t) exp(— icjt) dt 
which gives the power as |^. Here T is the duration of the observation. Because the signal has a small duty cy- 

cle instead of being purely sinusoidal, several harmonics are present in the Fourier power spectrum. One usually performs 
"harmonic summing" up to the -n}^ harmonic to increase the search sensitivity by a factor of up to y^n (see, e.g., Lorimer & 
Kramer 2005). From Equation ((3]), the power of the m"' harmonic in the Fourier spectrum can be written as 
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where we note that do contributes only as a constant phase factor and 
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is the observed angular spin frequency of the m* ^ harmonic. As mojp changes with t, the power gets distributed over adjacent 
Fourier bins leading to a loss in sensitivity. Nevertheless, the pulsar can be detected in the Fourier bin where the power is 
maximum. In the case of a pulsar moving with a constant velocity, the power of the m* harmonic will be peaked at the 
Fourier bin of angular frequency muip(l + ^) which is the conventional Doppler effect. 

Following ,]K91 let us define three efficiency factors, 7im, 72m and 73m. The factor 71m is the ratio of the height of the 
power of the m*'^ harmonic in the Fourier spectrum when the acceleration, jerk, and other higher order derivatives of the 
pulsar are non-zero, to the height when these terms are zero, i.e., the pulsar has a constant velocity. This can be written as 
follows 
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Figure 1. The orbit of a pulsar in a binary system, projected onto a plane containing the line of sight, i.e., perpendicular to the sky 
plane. The semi-major axis of the projected ellipse = Opsini where i is the orbital inclination angle and ap is the semi-major axis of 
the actual orbit of the pulsar. P is the location of the pulsar at t = and P is the location of the pulsar at any arbitrary time t. f is 
the true anomaly and ro is the longitude of the periastron. do is the distance between the pulsar and the Earth at t = and d is the 
distance between the pulsar and the Earth at t. 



where Ot, is a free parameter. The pulsar will be detected for such a value of which maximizes 71m with apparent frequency 
of the m"' harmonic as mtOp = mujp {1 + av/c). Clearly, the maximum value of = 1 is possible only when the pulsar 
moves with a constant velocity, giving = viq. 

The above efficiency factor 71™ essentially describes the sensitivity loss of a standard pulsar search. Acceleration searches 
attempt to improve upon this efficiency by accounting for the Doppler shifting of the p ulsar signal during the obs ervation. A 
review of some of the various techniques that have been developed so far can be found in lLorimer fc Kramer! (|2005l ). The most 
common type of acceleration search assumes that the line-of-sight acceleration during the observation is a constant value. This 
is a good approximation for binary periods that are significantly longer than the survey integration time, and we consider 
these searches using the efficiency factor 72m- Specifically, 7!^ is the ratio of the height of the power of the m*'^ harmonic in 
the Fourier spectrum when the jerk and other higher order derivatives of the pulsar are non-zero, to the height when these 
terms are zero, i.e., the pulsar has a constant acceleration. Expressing this idea mathematically, we have 
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where aa, ctv are free parameters. The use of an acceleration search algorithm leads to the detection of the pulsar for such 
a set of values of aa and a^ which maximizes 72m. Clearly, the maximum value of 72m = 1 is possible only when the pulsar 
moves with a constant acceleration, giving a„ = vio and aa ~ aio/2\. 

Due to the additional computational requirements, searches involving the unknown jerk term are less commonly carried 
out. To investigate their relative efficiency, however, we define 73m such that 73^ is the ratio of the height of the power of 
the m**^ harmonic in the Fourier spectrum when the derivative of the jerk and other higher order derivatives of the pulsar are 
non-zero, to the height when these terms are zero, i.e., the pulsar has a constant jerk. With this in mind, we may write 
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(12) 



where Uj, aa, and a^ are free parameters. Such a search process would detect a pulsar with a set of values of aj, aa and q„ 
which maximizes 73m. Clearly, the maximum value of 73^ = 1 is possible only when the pulsar moves with a constant jerk, 
giving Q„ = vio, aa = ciio/2! and aj = jio/3!. 

To compute the relevant line-of-sight terms in the above expressions, we need to consider the shape of the pulsar orbit. 
We denote the semi-major axis of the pulsar orbit to be Up and the semi-major axis of the ellipse projected perpendicular to 
the plane of the sky to be = ttpsini where i is the orbital inclination angle. This projected ellipse is shown in Fig. [T] where 
/ is the true anomaly and zu is the longitude of the periastron. The pulsar is located at P at f = and at P at an arbitrary 
time t. If r is the radius vector of the position of the pulsar at t with magnitude r, then the projection of r along the line of 
sight is given by 



ri 



r sin (/ + vj) 



(13) 
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Similarly, the projection of the magnitude ro of the radius vector of the pulsar at t = can be written as 

rio = ro sin (/o + -co) , 
where fo is the true anomaly at t = 0. We can write 

/ vidt^ {ri - no). 
Jo 

So, Equation ([7} becomes 

d = do + {n - no) . 
Using Equation (fTS)) . Equations (fTO)) . ([11} and ([12} become 



7lm(Qi.,r) = 



exp 



(n — no — Out) 



72m(aa, a„,T) = — I / exp 



— [ri — r;o — aat — 



and 



73m(Qj,aa,a!„,r) = — 



exp 



/ 3 2 \ 

— (n-r;o-Qjt - Qat - Qt,tj 



dt\ 



For the general case of an orbit of eccentricity e, we have 

r = (l — e^) (1 + e cos Z)""*^ 

so that 

ro = ttp (1 - e^) (l + ecos/o)"V 
With this definition. Equations (|13p and (|14p become 



ri — Qp (l — e ) (1 + e cos /) sin (/ + uj) 



and 



(14) 

(15) 

(16) 

(17) 
(18) 

(19) 

(20) 
(21) 
(22) 
(23) 



rjo = Hp (l — e ) (1 + ecos/o) sin (/o + rj7) 
Now, differentiating Equation (|13p and using Equation (|20p we get the expression for the velocity along the line of sight as 

27r a. 



vi =ri = — 



[cos (/ + vj) + e cos (zu) 



(24) 



where Po is the orbital period of the pulsar and the orbital angular frequency is ujo — 2tt/Po. The quantities Sj, and Po are 
related by Kepler's third law as 

1/3 



Mo 



a„ — ap sm i = an ■ 



^ ) G (Mp + M,) 



Mo 



Mp + Mo 



(25) 



Mp + Mo 

where Mp is the mass of the pulsar and Mo is the mass of the companion and G is Newton's gravitational constant. 
Differentiating Equation (|24|) . we get the expression for the line-of-sight acceleration 

2 ' 

■sin(/ + cc7) (1 + ecos/)^ (26) 



/2^\ 
a; = -Ui = - TT 



Po) (1 



o2^2 



Similarly, differentiating Equation (|26|) . we get the expression for the line-of-sight jerk 

3 ' 

^3 I 



Ji 



27r 



„2^7/2 



(1 -f e cos /) [cos (/ -I- vj) -f e cos [vj) — 3e sin (/ -I- ro) sin (/)] 



(27) 



In principle, one can continue differentiating to get higher order derivatives (e.g. the "jounce" is the time derivative of the 
jerk). The expressions for vi and ai were previously derived by Freire. Kramer, fc Lvnel ( 2001 . 2002). Assigning / = /q in 
Equations p4|l . (|26|l and (|27|l we get the values of vio, aio and jio respectively. 

To compute Equations (|17|1 . (I18|l and (I19p numerically, we need to solve Kepler's equations given below 

M = ujo{t - Tp) (28a) 



E - esinE ^ M 
f = 2tan"^ 



1 + £■ 
tan — 

1 - e 2 



(28b) 
(28c) 



where M is the mean anomaly, E is the eccentric anomaly and Tp is the epoch of the periastron passage. Mo — —uJoTp, Mo 
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Figure 2. Variation of the line-of-sight velocity with time over a complete orbit for a binary pulsar with Mp = 1.4 Mq, Mc = 0.3 Mq, 
Po = 0.5 day, i = 60°, Tp = 0, e = 0.5 for different ro, and the line-of-sight velocity for the same binary having zero eccentricity. 

being the mean anomaly at t = 0. For a circular otbit, e = 0, -ou — 0, f = E — M, giving / = ujot -f /o, /o = Mo — —uJoTp and 
we can get the expressions given by JK91 by substituting A = apUJo and (f> = fo + n/2 (or if = 270°, (f> = /o). 

In Fig.[2l we plot the line-of-sight velocity with respect to time over a complete orbit for a binary pulsar with Alp = 1.4 Mq, 
Mc = 0.3 A/0, Po = 0.5 day, i — 60°, Tp = 0, e — 0.5 for different zu between 0° — 360° and compare with the line-of-sight 
velocity plot in case of zero eccentricity. Like the circular case, line-of-sight velocity curve for the eccentric orbit is also 
symmetric over the half orbital period when zu is either 0°, 90°, 180° or 270° (upper panel), but not for any other values of 
zu (lower panel). In the lower panel of Fig. [21 is chosen to be 60°, 120°, 240° and 300°, each in a different quadrant, but 
have the same average line-of-sight velocity over the entire orbit. The same can be found for other values of zu. That is why, 
in the next section, when we perform our analysis for different zu, we choose values only in the ranges of 0° — 90°. The high 
value of eccentricity chosen here is justifiable as binary pulsars with eccentricity 0.5 or even higher exist in globular clusters, 
mainly as a result of stellar encounters (for a review, see Camilo & Rasio 2004). 

For relativistic binaries, Po, e, and zu change with time, and the rates of change of these parameters are classified as 
Post-Keplarian parameters. These changes need to be considered while performing pulsar timing analysis which involves pulse 
arrival times over a long span of time (days or even years); but, for the present study, we are concerned about the observability 
of pulsars, and these changes are negligible over the duration of any particular observation. 

3 ANALYSIS 

To demonstrate the impact of orbital motion on detectability, we will use the above framework to determine the maximum 
value of 7im, 72m and 73m. As different stellar and orbital parameters like Pp, Po, Mc, Mp, sini, zu come into the equations 
given in Section [51 we perform our analysis for different values of these parameters. For a particular binary, the above 
parameters are fixed as the change in these parameters can be neglected over the duration of the observation. The orbital 
phase, i.e. the true anomaly /, is, however, not known a priori, and changes during the observation. 

We divide the full-orbit (because of the asymmetry of an eccentric orbit which we have mentioned earlier) in equal steps 
of 10°, having fo.qS as 0°, 10°, 20° . . . 350 {q going from to 36), where /o,qS are the values of the true anomaly at the start of 
the observation. For each fo^q, we calculate -Eo,q, -Mo, 5 and Tp^q using Kepler's equations (|28al [28b1 128c|) . We use these Tp^qS 
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Table 1. Values of 7mli for different values of /o (true anomaly at the begining of the observation) and corresponding values of w for a 



binary with Mp = 1.4 Mq, Mc = 0.3 Mq, j = 60°, ro = 0° 



: 0.5, Pp = 0.01 s, m = 4, T = 500 s. 



/o 


71 m 


w 


0° 


0.28 


0.111 


10° 


0.23 


0.112 


40° 


0.19 


0.131 


50° 


0.19 


0.143 


60° 


0.21 


0.160 


70° 


0.22 


0.182 


170° 


0.98 


0.970 


180° 


0.99 


1.000 


190° 


0.86 


0.970 


350° 


0.38 


0.112 




Figure 3. Variation of the line-of-sight acceleration with true anomaly over a complete orbit for a binary pulsar with Mp = 1.4 Mq, 
Mc = 0.3 Mq, i = 60°, ro = 0°, e = 0.5, Po = 0.1 day. 



to calculate M, E and / at any time t (between and T) again solving equations H28al [28bl I28cp . This takes care of the 
change in orbital phase during the observation. Then we use Equations (|24p . H26p and (|27p to get — vi(t), aa ~ ai{t)/2\ 
and aj = ji{t)/3\ for different choices of t between and T and, for each set of a„, and aj, we perform the integrations 



given in Equations (|17p . (|18|) and (|19p numerically. The maximum values of of these integrations are the desired values 7mii3, 
7m2g and 7^39 respectiveljij. After this stage, we perform the weighted average of these efficiency factors over different fo.q. 

To perform this weighted average, we need to calculate the probability of each fo,q using the fact that the probability 
of a pulsar to be at a particular position in the orbit is directly proportional to the time it spends at that position. As the 
areal velocity is constant (Kepler's second law), the ratio of the time spent at true anomalies /i and /2 can be expressed as 
(1 + e cos : (1 + e cos /2)~^. So the weight factor Wq for each fo,q is given by = (1 + e cos fo,q)~^ / (1 + e cos /o.ig)"^. 

/o,i9 = 180° corresponds to the apastron, where the pulsar spends most of its time. Using the values of 7migS, ^m2qS and 
ImSijS for each fo,q, weighted averages can be estimated as: 



E; 'ym2,avq — ^ — ^ ; rrrLo.tLuu — — ^ 

q ^1 E9 ^1 Wq 

In this paper, we always report these average efficiency factors unless otherwise stated. We will henceforth skip the subscript 
avg, in section U and subsequent sections. 

In Table [T] we show the values of 7mi<j for different values of fo and corresponding values of w for a binary with 
Mp = 1.4 Mq, Mc = 0.3 Mq, i = 60°, uj = 0° , e = 0.5, Po = 0.1 day, Pp = 0.01 s, m = 4, T = 500 s. From Fig.[3 we see that 
ai has the largest magnitude at / = 57°. Both fo — 40° and fo = 50° contains this value of / during the observation (the 
pulsar moves 20.83° during a 500 s long observation), so we get the smallest value of ^mii here. Note that although fo — 170° 
and fo — 190° are equally likely (same value of w) and correspond to the same value of a;, but as the values of a; during 
next 500 s are different, they give different values for 71m ■ Similarly, both f — 0° and / = 180° correspond to a; — 0, but the 
values of ai during next 500 sec after fo — 0° and fo — 180° are different, so they result different values of 7im- 



J2q 7m2? 



Eg W;<7 "fmSq 



g= 1,2...36. 



(29) 



Remember that only comes in the expression for 7mig , and only aa and come in the expression of 7m2i3 . 
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4 RESULTS 

4.1 Neutron Star - White Dwarf (NS-WD) Binaries 

If the neutron star is a recycled milhsecond pulsar, then the orbit is ex pected to be circular because of the strong tidal phas e 
the past when the neutron star accreted matter from the companion (jAlpar et al ]| 19821 : iRadhakrishnan fc Srinivasanll 19821 ) . 



m 



However, in dense stellar environments like globular clusters or nuclear star clusters near the Galactic center, stellar interactions 
lead to eccentric orbits for recycled pulsar binaries. Many such systems in globular clusters are known at presen10. Moreover, 
three body interactions can also produce millisecond pulsars in eccentric orbit. One such system is PSR J1 903+0327, which 
was a member of a triple system i n the past, but has become a binary by ejecting one member of the system ij Champion et al 



2008 



3; Portegies Zwart et al. 2011). Fig. [2] shows the line-of-sight velocity plots for a canonical example of this category. 



It is clear from Equation ((9]), that the spread of the power among adjacent Fourier bins is larger for higher values of m 
or larger values of iUp, thus the values of the efficiency factors decrease for higher m or smaller Pp. Similarly, a smaller value 
of Po gives higher values of vi , a; and ji and so on (all other parameters being fixed) resulting in a larger spread of the power 
among adjacent Fourier bins, i.e., reduction of the efficiency factors. It is evident from Eqn. (|25p that higher value of sini will 
increase the values of vi, ai, ji and so on, leading to a decrease in the efficiency factors. To demonstrate these facts graphically, 
we show variations of 71m with with Po and Ps for different parameters. In all the subsequent plots, the X-axis shows the 
values of Pp, the Y-axis shows the values of Po and the color code represents the values of 71m. The contour for 71m — 0.5 
is also shown in each plot. We can say that the pulsars having 71m < 0.5 (lying on the bottom-left side of the 71m ~ 0.5 
contour) are difficult to detect. Fig. |4] shows the variation for different values of zu. We see that the 71m = 0.5 contour shifts 
rightward with the increase of w, making a larger portion of the phase-space difficult to detect. Fig. [5] shows the variation for 
different values of i, and as expected, the increase of i makes a larger portion of the phase-space difficult to detect. In Fig. [6l 
we show variations of 7i,„ with with Po and Ps for different Mc keeping all other parameters fixed, and we see that the values 
of 7im are lower for higher values of Mc when all other parameters are the same. This happens because Op increases with 
the increase of Mc (if Po remains the same) giving higher values of vi, ai and ji etc. resulting in a larger spread of the power 
among adjacent Fourier bins, i.e., reduction of the efficiency factors. As an example, if Ale increases from 0.1 Mq to 0.8 Mq, 
the fractional increase in Op is 5.2. On the other hand, from Fig [71 we have not seen any significant effect of the variation of 
Mp in the realistic range of 1 — 2 Mq, because the fractional increase of Up is only 0.32 when Mp increases from 1 Mq to 
2 Mq . In Fig. [SI we show variations of 71m with with Po and Ps for different values of e and find that it does not play any 
significant role unless very high when 71m increases with the increase of e. In Fig. [9] we show the variations of 71m with with 
Po and Ps for m = 1 and m = 7 keeping all the other parameters fixed and, as expected, 71m decreases for higher values 
of m. In Fig. IIOI we show variations of 71m with with Po and Ps for different values of T and find that a smaller value of 



T mak es a larger portion of the phase space easy to detect. This fact was also noticed bv iRamachandran fc Portegies Zwart 



1 1998t ). On the othe r hand, smaller T reduce s the value of minimum detectable fiux density for a pulsar, independent of its 



orbital parameters (|Lorimer fc Kramerll2005l l. Hence an optimal strategy is needed to choose a value of T depending upon 
the motivation of the survey and efficiency of the telescope. 

In Fig. 111! we show variation of 72m and 73m with Po and Ps keeping = 30°, i — 60° Alp — 1.4 Mq, Ale = 0.3 Mq, 



m = 4 and T = 1000 s. These can be compared to Fig. 4(b) which shows the variation of 71m for the same set of parameters. 
The increases in efficiency factors by using acceleration or acceleration-jerk searches are significant. We notice that the slope 
of the contour is different and the use of acceleration or acceleration-jerk searches are more effective for short orbital periods. 
This result agrees with that obtained by JK91. 

In Table [3l we present the efficiency factors for T = 1000 s and m = 4, averaged over different /qs as described earlier for 
the NS-WD binaries for which all the required parameters, e.g., Pp, e, Po, sini, -oo, Mc and Mp are known, for few cases, the 
values of Alp are not known, so assumed to be 1.4 AIq. Table [2] shows the parameters for these binaries and corresponding 
references. 

For the sake of simplicity, in most cases, we compute the efficiency factors for the m = 4 harmonic, but the procedure 
will be the same for any other harmonic. The value of the efficiency factor at any harmonic m2 for spin period Pp^ would be 
the same as that at the harmonic mi and spin period Pp^ if Pp2 ~ ^Ppi if we keep all other parameters unchanged. This 
fact can be used to read the efficiency factors for any other harmonic from the plots. 

In the case of a r eal pulsar, the number of harmonics present in the Fourier spectrum of a pulse is roughly reciprocal of 
the pulse duty cycle ( Lorimer fc Kramer 2005l l . Harmonic summing is used to get the power from each of these harmonics. 



The power of lower harmonics are higher as well as the efficiency factors. As a combined effect, the contributions from 
higher harmonics bec ome more significant for a binary pulsar. An an example, we fit the pulse profile of PSR J1802— 2124 



( Faulkner et al\\2004 ) with a Gaussian, and compute the Fourier spectrum. Table [4] shows the values of 71m for first 10 



harmonics, power {Dm) of each harmonic (in arbitrary unit) and degraded power (71m Dm) of each harmonic for different 
^ http: / / www.naic.edu / '^pfreire/GCpsr.html 
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(c) ro = 60° (d) -w = 90° 

Figure 4. Variation of witii Po and Ps for different ro. For eacli case, i = 60°, Mp = 1.4 Mq, Mc = 0.3 Mq, e = 0.5, m = 4, 
T = 1000 s. 

observation duration, e.g. T =500, 100, 1500 and 2000 s. As expected, for any fixed value of T, 71^ decreases with increasing 
value of m, and for any fixed value of m, 71^ decreases with the increase of T. Note that, for any fixed value of T, the decrease 
in the value of 7i„ Dm from m = 1 to m = 10 is slower than that of Dm ■ 



4.2 Double Neutron Star (DNS) Binaries 

Presently there are eight confirmed double neutron star (DNS) binaries and four more candidate DNSs. Among these twelve, 
all required parameters to calculate efficiency factors are known only for four. We present the efficiency factors for these DNSs 
in Table [6l for T — 1000 s, m = 4. The parameters for these DNSs with corresponding references are given in Table [5] 

In Fig. [12] we show the line-of-sight velocity curves for a hypothetical DNS having Adp = 1.35 Mq, Mc — 1.25 Mq, 
Po = 0.5 day, i = 60°, Tp = 0, e — 0.5 for different zu. DNS systems can have such high eccentricity and short orbital period. 
One example is the Hulse- Taylor binary PSR B1913+16, which has e — 0.617 and Po = 0.323 day. Comparing Fig. 1121 with 
Fig. [2I we see that the amplitudes of radial velocity curves are ~ 3 times larger for a DNS in comparison with a NS-WD binary 
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(a) i = 30° (b) i = 90° 



Figure 5. Variation of 71m with Po and Ps for different i. For each case, = 60°, Mp = 1.4 Mq, Mc = 0.3 Mq, e = 0.5, m = 4, 
T = 1000 s. 



Table 2. Parameters for NS-WD binaries for which we can calculate efficiency factors. The columns from left to right show the name of 
the pulsar, spin period, orbital period, orbital eccentricity, sine of the inclination angle, longitude of periastron, companion mass, pulsar 
mass and corresponding references. 



Pulsar 




Po 


e 


sin i 




Mc 


Mp 


Refs. 




(sec) 


(day) 






(dog) 


{Mq) 


(.Mq) 




J0437-4715 


0.005757 


5.74105 


1.918 xlO"^ 


0.674 


1.222 


0.254 


1.76 


Verbiest et al. f2008~l 


J0751+1807 


0.003479 


0.26314 


5.0 xlO"'^ 


0.9121 


45.0 


0.191 


1.26 


Nice et al. ("2005); ATNF catlogue* for 


















Nice, Stairs, & Kasian (2007. for M„) 


J1600-3053 


0.003598 


14.34846 


17.369 xlO-5 


0.8 


181.85 


0.6 


1.4* 


Verbiest et al. (2009) 


J1603-7202 


0.014842 


6.30863 


9.3 xlO"'' 


0.89 


169 


0.14 


1.4* 


Hotan, Bailes. & Ord (2006) 


J1614-2230 


0.003151 


8.68662 


1.30 xlO"*^ 


0.9999 


175.1 


0.500 


1.97 


Demorest et al. (2010) 


J1640+2224 


0.003163 


175.46066 


7.973 xlO~* 


0.99 


50.731 


0.15 


1.4* 


Lohmer et al. (2005) 


J1713+0747 


0.00457 


67.82513 


7.494 xlO"^ 


0.95 


176.191 


0.28 


1.3 


Splaver et al. (2005) 


J1738+0333 


0.00585 


0.35479 


3.4 xlO"'^ 


0.5388 


155.695 


0.19 


1.46 


Freire et al. (2012) 


J1802-2124 


0.012648 


0.69889 


2.48 xlO""^ 


0.9845 


20.0 


0.78 


1.24 


Ferdman et al. (2010) 


J1903+0327 


0.00215 


95.17412 


0.4367 


0.976 


141.652 


1.03 


1.67 


Freire et al. (2011) 



t : |http://www. atnf.csiro.au/research/pulsar/psrcat/expert. html| 
* : Mp not measured, taken as 1.4 Mq 



Table 3. Efficiency factors and corresponding parameters for NS-WD binaries (for which all required parameters are known) for T = 1000 
s, m = 4 (averaged over different TpS as described in the text). 



pulsar 


71 m 


72m 


73m 


J0437-4715 


1.00 


1.00 


1.00 


J0751+1807 


0.29 


0.64 


0.98 


J1600-3053 


1.00 


1.00 


1.00 


J1603-7202 


1.00 


1.00 


1.00 


J1614-2230 


0.99 


1.00 


1.00 


J1640+2224 


1.00 


1.00 


1.00 


J1713+0747 


1.00 


1.00 


1.00 


J1738+0333 


0.56 


0.95 


1.00 


J1802-2124 


0.50 


0.98 


1.00 


J1903+0327 


1.00 


1.00 


1.00 
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0.6 



0.2 




(a) A/c = 0.01 Mq 



(b) Mc = 0.1 Mq 




0.8 



0.6 



Pp(s) 
(c) Mc = 0.8 M0 

Figure 6. Variation of 71^ with Po and Pj, for different Mc. For eacli case, i = 60°, = 60°, Afp = 1.4 Af©, e = 0.5, m = 4, T = 1000 s. 

Table 4. Values of 71^ for first 10 harmonics of PSR J1802— 2124, power (Dm) of each harmonic (in arbitrary unit) and degraded power 
(7im ^"i) °f ^^^^ harmonic for T =500, 100, 1500 and 2000 s. The pulse profile was fitted with a Gaussian. 



harmonic number (m) power (Dm) 7in 



= 500 

7?„ 



T = 1000 s 

7lm 7?m -Dr, 



T = 1500 s 



■yir. 



7i„ 



7lr 



2000 s 
7L ^" 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



1.00000 
0.88965 
0.62585 
0.34728 
0.15156 
0.05199 
0.01412 
0.00310 
0.00056 
0.00008 



0.99 
0.97 
0.93 
0.87 
0.81 
0.74 
0.70 
0.67 
0.64 
0.61 



0.982998 
0.830733 
0.536621 
0.264492 
0.099354 
0.028554 
0.006924 
0.001394 
0.000230 
0.000029 



0.87 
0.67 
0.57 
0.50 
0.46 
0.42 
0.39 
0.38 
0.35 
0.34 



0.761651 
0.400844 
0.203214 
0.088176 
0.031453 
0.009198 
0.002202 
0.000441 
0.000070 
0.000009 



0.64 
0.48 
0.41 
0.36 
0.32 
0.30 
0.27 
0.25 
0.24 
0.23 



0.412938 
0.206042 
0.103509 
0.043970 
0.015686 
0.004363 
0.001009 
0.000196 
0.000032 
0.000004 



0.51 
0.38 
0.31 
0.27 
0.25 
0.23 
0.22 
0.20 
0.19 
0.18 



0.257970 
0.128296 
0.060903 
0.025762 
0.009186 
0.002693 
0.000674 
0.000125 
0.000020 
0.000003 
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(a) Mp = I.OMq (b) Mp = 2.0 Mq 

Figure 7. Variation of 71^ with Pq and Ps for different Mp. For eacli case, i = 60°, ro = 60°, = 0.3 Mq, e = 0.5, m = 4, T = 1000 s. 



Table 5. Parameters for tlie DNSs for which we can calculate efficiency factors. The columns from left to right show the name of the 
pulsar, spin period, orbital period, orbital eccentricity, sine of the inclination angle, longitude of periastron, companion mass, pulsar mass 
and corresponding references. 



DNS 




Po 


e 


sini 




Mc 


Mp 


Refs. 




(sec) 


(day) 






(deg) 


(Mq) 


{Mq) 




J0737-3039A 


0.022699 


0.10225 


0.08778 


0.9997 


87.033 


1.249(B) 


1.338(A) 


Kramer et al. C2006) 


J0737-3039B 


2.773461 


0.10225 


0.08778 


0.9997 


267.033 


1.338(A) 


1.249(B) 


Kramer et al. (2006) 


B1534+12 


0.037904 


0.42074 


0.27368 


0.975 


274.769 


1.35 


1.328 


Konacki et al. (2003); Stairs et al. (2002. for Mtnt) 


J1756-2251 


0.028461 


0.31963 


0.18057 


0.95 


327.825 


1.312 


1.258 


Ferdman (20081 


J1807-2500Bt 


0.004186 


9.95667 


0.74703 


0.996 


11.335 


1.206 


1.365 


Lynch et al. (2012) 


B1913+16 


0.059030 


0.32300 


0.61713 


0.71 


292.545 


1.389 


1.440 


Weisbere. Nice. & Taylor (2010"); 



^Damour_&_^hafer (1988, for sin i) 



t : candidate DNS. 



having the same values of Po, e, sini and vj. This fact leads to a reduction in the efficiency factors for DNSs in comparison 
with NS-WD binaries. 

In Fig 1131 we show the variations of 7im, 72m, and 73^ with Po and Pg for Mp — 1.35 Mq, M^ = 1.25 Mq, zu = 30°, 
i = 60°, e = 0.5, m = 4, and T = 1000 s. The variations of 71^ with Po and Ps for e = 0.1 and e = 0.8 keeping all other 
parameters unchanged are shown in Fig [14] Again, it is clear that a higher value of eccentricity is favou rable for detection, i.e. 
results a higher value of 71m. This fact was also observed by Ramachandran fc Portegies Zwart 1 19981 ). who studied eccentric 
systems only for DNS binaries. 



Table 6. Efficiency factors for DNSs (for which all required parameters are known) for T = 1000 s, m = 4. 



DNS 


71 m 


72m 


73m 


J0737-3039A 


0.21 


0.41 


0.56 


J0737-3039B 


0.98 


1.00 


1.00 


B1534+12 


0.52 


0.92 


1.00 


J1756-2251 


0.41 


0.85 


0.99 


J1807-2500B 


0.92 


1.00 


1.00 


B1913+16 


0.68 


0.92 


0.94 
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Figure 8. Variation of 7i„i with Po and Ps for different e. For eacfi case, i = 60°, ro = 60°, Mp = 1.4 Mq, Mc = 0.3 Mq, m = 4and 
T = 1000 s. 
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4.3 Neutron Star - Black Hole (NS-BH) Binaries 



Although no neutron star - (stellar mass) black hol e (NS-BH) binary is known at present, there is high possibility of detecting 
such systems with forthcoming facili ties like SKA ( Smits_et_aZj 200S ). 

A number of theoretical studies ( Pfahl. Podsiadlowski. fc Rappaport 20051 : Kiel &: Hurlevl[2C)09l: Faucher-Giguere &: Loebl 



2011 ) exist which predict probable parameters for these systems. The population synthesis model of lPfahl. Podsiadlowski. fc Rappaport 
l|2005l l predicts mi ldly recycled pulsars in short orbits (Po in the range of 0.004 — 0.8 days) with small eccentricities, 
while the study by Kiel fc Hurley predicts recycled pulsars with orbital periods ranging from 0.1 days to 10^ days. 

Faucher-Giguere fc Loebl (|2011 ) found that recycled pulsars with black hole companions in highly eccentric orbits with Po 



between 0.003 and 0.6 days can form due to stellar encounters in dense stellar environments. In this work, we calculate the 
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0.6 



0.2 



0.001 0.01 



(a) 72m (b) 73m 

Figure 11. Variations of 72m and 73m with Po and Ps keeping ro = 30°, i = 60° Mp = 1.4 Mq, Mc = 0.3 Mq, m = 4 and T = 1000 s. 
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Figure 12. Variation of the line-of-sight velocity with time over a complete orbit for a DNS, having Mp = 1.35 Mq, Mc = 1.25 Mq, 
Po = 0.5 day, i = 60°, Tp = 0, e = 0.5 for different zu, and the line-of-sight velocity for the same binary having zero eccentricity. 

efRciency factors for 0.1 < Po < 100 days and 0.001 < Pp < 1, as we chose for the cases of NS- WD and NS-NS binaries. We 
take the mass of the black hole as 10 Mq, which agrees with the observations (jFarr et aZ.ll201in . 



In Fig. 1151 we show the variation of the radial velocity with time for a NS-BH binary with Mp — 1.4 Mq, Mc = 10.0 Mq, 
Po = 0.5 day, i — 60°, Tp = 0, e = 0.5 for different zu, as well as the circular case (keeping all other parameters fixed). We see 
that the amplitudes of the line-of-sight velocity curves are ~ 3 times larger for a NS-BH in comparison with a NS-NS binary 
having the same values of Po, e, sini and w (or ~ 9 times larger for a NS-BH in comparison with a NS-WD binary). So it is 
clear that T — 1000 s is too high for this case and will lead to very small efficiency factors and we decided to use T — 500 s. 

In Fig. 1161 we show the variations of 71,71, 72m, and 73m with Po and Ps for vj = 30°, i = 60°, e = 0.5, m = 4, and 
T — 1000 s. The variations of 71m with Po and Ps for e = 0.1 and e = 0.8 keeping all other parameters unchanged are shown 
in Fig 1171 A higher value of eccentricity is favourable for detection, i.e. results a higher value of 71m • 



5 CONCLUSIONS 

We have generalized the earlier work of JK91 to calculate the signal-to-noise degradation of pulsars in binary orbits with 
arbitrary eccentricity. We have applied the framework to compute degradation factors for a variety of orbital configura- 
tions, and show that it can be quite substantial. We have also demonstrated how the degradation can be recovered by using 
acceleration search or acceleration-jerk search algorithms. The analysis should prove invaluable to a wide variety of popula- 
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tion studies since this work provides, for the first time, accurate accounting for this important effect. Software to calculate 
the degradation factors for arbitarary orbital parameters, harmonic summing and survey integrations is available online at 
|http://psrpop. phys.wvu.edu/binary, . 
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(a) e = 0.1 



(b) e = O.i 



Figure 14. Variation of 71^ with Po and Ps for Mp = 1.35 Mq, Mc = 1.25 Mq, to = 30° , i = 60° , m = 4, and T = 1000 s, for e = 0.1 
and e = 0.8. 
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Figure 15. Variation of the line-of-sight velocity with time over a complete orbit for a binary pulsar with a black hole companion with 
Mp = 1.4 Mq, Mc = 10.0 Mq, Po = 0.5 day, i = 60°, Tp = and e = 0.5 for different ro, and the line-of-sight velocity for the same 
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Figure 16. Variation of 71m, 72m, and 73m with Po and Ps for Mp = 1.4 Mq, Mc = 10.0 M©, ro = 30°, i = 60°, e = 0.5, m = 4, and 
T = 1000 s. 
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(a) e = 0.1 (b) e = 0.8 

Figure 17. Variation of 71™ with Po and Ps for Mp = 1.4 Mq, Ma = 10.0 Mq, m = 30°, i = 60°, m = 4, and T = 1000 s, for e = 0.1 
and e = 0.8. 
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